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Using first-principles calculations, we investigate the electronic structures and binding properties 
of nicotine and caffeine adsorbed on single- walled carbon nanotubes (CNTs). The purpose of this 
study is to determine whether CNTs are appropriate for sensing nicotine and caffeine molecules. 
We find that the two types of molecules exhibit very different binding characteristics on CNTs; 
nicotine binds to the CNT with an adsorption energy of 0.46 eV, whereas caffeine adsorbs with a 
higher binding energy of 0.88 eV. These differences are discussed by analyzing the modification of 
the electronic structure of the molecular-adsorbed CNTs and charge transfer with the adsorbates. 
Especially, we find a special binding feature between the caffeine molecule and the CNT formed by 
electron donation and back-donation. 



I. INTRODUCTION 

Carbon nanotubes (CNTs) have attracted much aca- 
demic and industrial interest because of their remarkable 
physical and chemical properties ™JJ One of the most 
intriguing characteristics of CNTs is a very large surface 
ares ji2jti4j comparable to that of activated carbon. This 
property has been exploited to develop various applica- 
tions including hydrogen storage, filters removing toxic 
compounds, and chemical sensors. However, CNTs were 
not regarded as a good storage medium for hydrogen due 
to their inert surfaces, where H2 molecules adsorb very 
weakly via van der Waals interactions^^ in spite of 
overestimated earlier studies {SMS The binding proper- 
ties of hazardous molecules on CNTs were also examined 
theoretically and experimentally!**"*^ One of the most 
promising applications of CNTs is a sensor application. 
CNT-based sensors were developed and demonstrated 
to detect a wide range of molecules such as ammonia 
(NH^EZI nitrogen dioxide (NO 2 )| P 3 oxygen (0 2 )!^alco- 
holpSand other molecules JSHUl They verified that CNT- 
based sensors are able to detect not only the type of gas 
molecule, but also the concentration of the adsorbate. 

Nicotine and caffeine are both alkaloids and therapeu- 
tic compounds. They serve as stimulants for the central 
nervous system in humans. Through tobacco (cigarette) 
and coffee, nicotine and caffeine are the two addictive 
drugs used heavily by humans. Because of their psy- 
choactive effects, nicotine and caffeine could have posi- 
tive effects. For example, nicitone appears to enhance 
concentration and memory due to the increase of acetyl- 
choline,^ and caf feine seems to enhance performance in 
endurance sports ! 34 * 35 ! However, these drugs often pro- 
duce side effects. Both of them may increase the heart 
rate that further limits the bodys ability to maintain 
homeostasis during exercise. Overdose of these drugs is 
severe. Especially, nicotine overdose can be deadly, and 
excessive ingestion of caffeine over extended periods of 
time results in a toxic condition (caffeinism), with symp- 
toms of vomiting, elevated blood pressure, rapid breath- 
ing, heart palpitations, and insomnia. Therefore, the fab- 



rication of sensors and filters with high sorption capabil- 
ity of these chemicals is desired for public health. 

In this paper, we report our study on the adsorption 
properties of nicotine and caffeine on CNTs. We ex- 
amined the modification of the electronic structures of 
CNTs due to adsorbates. We found that both nicotine 
and caffeine have non-covalent interaction with the CNT 
through 7r-stacking. Although nicotine and caffeine are 
similar chemical compounds, containing aromatic rings 
and nitrogen atoms, it was found that they show very 
different binding geometries and binding energies. 



II. COMPUTATIONAL DETAILS 

We performed first-pri nciple s calculations based on the 
density functional theory^EHto study the structural and 
electronic properties of the adsorbates on CNTs. We 
employed a plane- wave basis set and the projector aug- 
mented wave (PAW) implemente d in the Vienna ab ini- 
tio Simulation Package (VASP)P^ The cutoff energy 
for the kinetic energy was set to 500 eV. The generalized 
gradient approximation (GGA^PHl anc [ the localized den- 
sity approximation (LDA)p^ were used to describe the 
exchange-correlation (XC) energy functional. In general, 
chemical bonding and electron transfer are well described 
within the LDA or the GGA. However, the dispersion 
or van der Waals (vdW) forces are not represented well 
by these two XC functionals because quantum electronic 
interactions in the regions of low electron densities are 
not correctly expressed. In these calculations, we con- 
sidered the correction of the vdW interactions between 
the adsorbates and the CNT, using Grimme's method,^ 
only for GGA. We also calculated the same model struc- 
tures with the LDA scheme, and compared its results 
with the those based on the GGA plus vdW. Since the 
LDA and the GGA plus vdW show the similar electronic 
structures, we present the only results obtained from the 
GGA plus vdW. 

To verify the validity of our computaional results, 
we also repeated our work with the linear combination 
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FIG. 1: Model structures of (a) nicotine and (b) caffeine. 



of pseudoato mic o rbitals method implemented in the 
SIESTA codeSHHU We used 210 Ry as the mesh cutoff 
energy and double-^ basis with polarization for the ba- 
sis set. All the model structures were relaxed until none 
of the residual Hellmann-Feynman forces acting on any 
atom would exceed 0.03 eV/ae, where <2b is the Bohr ra- 
dius. The binding energies were corrected using the basis 
set superposition error correction with ghost atoms. The 
results from both packages were almost identical and con- 
firm that our results are valid with DFT. 

We selected a zigzag (8,0) CNT, which is a semicon- 
ducting CNT, as the host material for molecular adsorp- 
tion. We used a tetragonal supercell with a length of 
16.94 A along the CNT axis and 20 A along the other 
two directions. We utilized 1x1x4 Monkhorst-Pack 
sampling 45 for the Brillouin zone integrations. On the 
other hand, a cubic supercell with a length of 20 A along 
all three directions and only 1 k point at T point were 
used for molecular calculations. Charge transfer between 
the adsorbate and the host CNT was estimated using the 
Mulliken population analysis. 



III. RESULTS AND DISCUSSION 

To begin with, we calculated the individual nicotine 
and caffeine molecules in vacuum condition to optimize 
the gas phase geometry. Figures [IJa) and (b) show the 
optimized geometries of nicotine and caffeine molecules. 
All the nitrogen atoms form essentially a planar structure 
with carbon atoms implying sp 2 orbital hybridizations, 
maintaining the aromatic characteristics. Figures [2ja)- 
(c) show the density of states (DOS), the highest occu- 
pied molecular orbital (HOMO) and the lowest unoccu- 
pied molecular orbital (LUMO) of the individual nicotine 
molecule. Similarly, figures [2] (d)-(f) represent the same 
quantities of the caffeine molecule. The energy gaps be- 
tween the HOMO and the LUMO levels are 3.34 eV and 
3.37 eV for the nicotine and caffeine molecules, respec- 
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FIG. 2: (a) DOS, (b) HOMO and (c) LUMO orbitals of nico- 
tine, (d) DOS, (e) HOMO and (f) LUMO orbitals of caffeine. 



tively, which are very similar to each other. However, 
the energy difference between the LUMO and LUMO+1 
(the state one level above the LUMO) of the nicotine 
molecule shown in figure |2|a) is much smaller that that 
of the caffeine molecule seen in figure |2jd). Interestingly, 
the HOMO wavefunction of the nicotine molecule is dom- 
inant at the five-membered ring as shown in figure [2jb), 
whereas that for the LUMO is placed at the phenyl group 
(six-membered ring) as displayed in figure |2jc). In con- 
trast, the HOMO and LUMO wavefunctions of the caf- 
feine molecule appear to be distributed over the whole 
molecule as shown in figures [2|e) and (f). 

Figure [3] shows the optimized structures of the (8,0) 
CNT with an adsorbed nicotine and a caffeine molecules. 
As shown in figure |3|a), the nicotine molecule prefers 
to contact with the CNT through its six-membered ring 
rather than its five-memebered ring. Moreover, the 
most stable binding gemoetry appears similar to Bernal 
("AB") stacking of graphite with the binding energy of 
Eh & 0.46 eV and the binding distance of ~ 3.29 A. 
On the other hand, the caffeine molecule is bound to 
the CNT stronger than the nicotine molecule by almost 
twice. Its binding energy and distance are « 0.88 eV 
and db ~ 3.25 A, respectively. In contrast to AB stacking 
in the nicotine adsorption, the caffeine molecule prefers 
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FIG. 3: Optimized structure of the (8,0) CNT with (a) 
an adsorbed nicotine molecule and (b) an adsorbed caffeine 
molecule. The figures on the left hand side show side view, 
whereas the right hand side figures do top view. To avoid 
confusion, carbon atoms in the CNT are shown yellow. 



"AA" to AB stacking, which is only ~ 5 meV less stable 
though. 

Figure [4] displays the electronic structures of the bare, 
nicotine-, and caffeine-adsorbed (8,0) CNTs. The bare 
(8,0) CNT has an energy band gap of 0.58 eV (fig- 
ure Qb)), whereas both nicotine- and caffeine- adsorbed 
CNTs increase their band gap to ~ 0.61 eV, as seen in 
figures [4^ a) and (c). In figures |4^a)-(c), we did not find 
a significant shift in the Fermi level by the molecule ad- 
sorption implying no significant charge transfer between 
the adsorbate and the CNT. In each of the band structure 
of the molecule-adsorbed CNT, there are two flat bands 
deemed to come from the molecule. For the nicotine 
(caffeine) case, such fiat bands are located near 0.6 eV 
(0.6 eV) and 1.35 eV (0.9 eV) below the Fermi level. Fig- 
ures [4^d) and (e) display the wavefunctions correspond- 
ing to the states labeld as "A" and "B" in figure |4ja) 
at the T point. The state A in figure Qd) shows a lo- 
calized state corresponding to the HOMO of the nicotine 
molecule, whereas the state B exhibits a small delocalized 
7r state on the CNT together with the nicotine HOMO- 
1 state (the state just below the HOMO) resulting in 
weak hybridization. Similarly, we show in figures [4jf) 
and (g) the wavefunctions plotted at the F point for the 
two states labeled as "C" and "D" in figure [4jc). Inter- 
estingly, due to the stronger interaction of the caffeine 
molecule with the CNT than that of the nicotine, the 
states C and D corresponding to the HOMO and the 
HOMO-1 states of the caffeine molecule show an interac- 
tion between an extended state of the CNT and a molec- 
ular orbital of caffeine. 

For further analysis, we calculated charge density dif- 
ferences shown in figure [5|a). Iso-density surfaces indi- 
cate regions of electron accumulation (red) and deple- 
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FIG. 4: The band structures of (a) the nicotine-adsorbed 
CNT, (b) the bare (8, 0) CNT and (c) the caffeine-adsorbed 
CNT. (d-e) The states A and B indicate wavefunctions at la- 
beled "A" and "B" in (a), (f-g) The states C and D indicate 
wavefunctions at labeled "C" and "D" in (c). The red and 
blue surfaces indicate the phase of the wave functions. 

tion (blue). From these plots, we also confirmed that 
there is no significant charge transfer between the ad- 
sorbate and the CNT. But it is shown that the nicotine 
molecule creates a local polarized regon at the CNT near 
the adsorbate as shown in figure [5|a), We found a more 
interesting feature in the caffeine-adsorbed CNT, where 
the caffeine molecule generates an electron depletion re- 
gion near the one end of the adsorbate, and at the same 
time an electron accumulation region at the other end of 
the adsorbate. This phenomenon can be interpreted as 
a "donation and back-donation" between the CNT and 
the caffeine molecule. 



IV. CONCLUSION 

In conclusion, we investigated the electronic struc- 
tures and binding properties of nicotine and caffeine 
molecules on (8,8) single-walled carbon nanotubes us- 
ing first-principles calculations. We found that both ad- 
sorbates show very different binding characteristics on 
CNTs. The nicotine molecule binds to the CNT with 
an adsorption energy of 0.46 eV, whereas the caffeine 
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molecule adsorbs with a higher binding energy of 0.88 
eV. To figure out the difference, we discussed the modi- 
fication of the electronic structures of the nicotine- and 
caffeine-adsorbed CNTs. Although there is no signifi- 
cant charge transfer between the CNT and the adsorbate, 
we found that there is an noticeable interaction between 
them. In particular, the caffeine-adsorbed CNT exhibits 
a special binding feature formed by electron donation and 
back-donation. 
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